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Abstract 

We perform molecular dynamics simulations of friction for atomically thin 
Xe films sliding on Ag(lll). We determine the inverse of the coefficient of 
friction (i.e slip time) by direct calculation of the decay of the center of mass 
velocity after applying an external force, as well as from the velocity auto- 
correlation function. We find that the slip time exhibits a drop followed by 
a sharp increase in a range of coverage near one monolayer. The slip time 
then levels off with further coverage increases in agreement with previously 
reported experiments. Our simulations suggest that the friction found in this 

system is dominated by phonon excitations. 
PACS numbers: 79.20.Rf, 64.60.Ht, 68.35. Rh 
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Although tribology, the study of friction and wear, has been of technological interest 
since ancient times the topic continues to rouse interest today P-p!2|. Rapid progress in 
experimental, theoretical and computational methods provides new insights into the atomic 
origins of frictional energy dissipation. When a thin film slides on a metal substrate there 
exists dissipation of energy via two mechanisms: (i) electronic excitations in the metallic 
substrate , and (ii) phonon excitations in the film or in the substrate [|l^ . The dissipation 



of energy can be characterized by the slip time r (i.e. the time it takes for the film's speed 
to fall to 1/e of its original value, assuming it is stopped by friction) or equivalently by a 
damping coefficient ~ 1/r. 

In this letter we study the phonon contribution to friction for Xe films sliding along a 
Ag(lll) substrate using molecular dynamics simulations. It is of great interest to find a way 
to determine the relative contributions of the phonon and electron contributions to friction, 
since to date it is not clear which is dominant. To this end, we compare our results with 
the slip-time versus coverage data reported by Daly and Krim 0, and with the electrical 
resistivity versus coverage data of Dayo and Krim 0]. The results of this comparison suggest 
that phonon excitations make a dominant contribution to the friction. 

We determine the slip time as function of coverage, defined as the number of atoms in 
the film per unit area. We treat a range of film coverages, from submonolayer to bilayer. 
The slip time is determined by two methods. In the first method, an initial center of mass 
velocity Vq is produced by an external force exerted on the film for t < 0. The external force 
is turned off for t > 0, and r is determined by the resulting velocity decay Voe^*/"^ . In the 
second method, no external force is applied. The slip time is determined by the behavior 
of the thermal equilibrium autocorrelation for the center of mass velocity as a function 
of time. The autocorrelation function is related to the linear response of the system to 
a small perturbation by Kubo theory [|13|. Linear response theory holds in the limit 
of zero driving velocity and zero response force. Therefore, this method is advantageous 
for describing experiments using the Quartz Crystal Microbalance (QCM ) |^,|^, since the 
QCM drives the film out of equilibrium only to a very small extent. This new method for 



calculating the slip-time is more appropriate for the small sliding velocities that occur in 
the actual experiment than previously used methods. Previous methods apply unphysically 
large forces (or, alternatively, shake the substrate at unphysically high frequencies) in order 
to generate velocities which are several orders of magnitude larger than those which occurred 
in the experiment. 

Molecular dynamics simulations of sliding friction have traditionally employed ther- 
mostats which add or remove energy to or from the system in order to achieve constant 
temperature. Thermostats used in molecular dynamics have the side effect of damping the 



atomic motions, resulting in an additional contribution to friction ||TT],[TO|. In order to study 
friction in a situation free of such complications, we employ a thermostat solely to establish 
statistical thermal equilibrium, but then turn it off while monitoring the system's properties. 
The absence of a thermostat in the simulation, allows us to focus exclusively on the film's 
dissipation due to phonons. 

The model Hamiltonian for film atoms of mass m at positions [k = 1, .., A^) is given 

by 

N 2 

^-E|^ + ^K---,r^), (1) 

k=l 

where is the momentum of the atom k, and the total potential U{ri, ■ ■ ■ , r^r) is given by 

N N 

f/(ri,---,rjv) = E^^(rfc)+ E ^(k, " r^l). (2) 

fc=l j<k=l 

Here, f/s(rfc) is a single particle potential describing the interaction between the k-th film 
atom and the substrate, and V{\rj — rk\) is the pair potential interaction between the j-th 
and k-th atoms in the film. 

The interaction between two Xe atoms is given by a Lennard- Jones potential 



V{r) = 4e 



(3) 



where e = 19.83 meV^, and a = 4.055 A |T^. The interaction between a Xe atom and the 
substrate can be described by JlSl a substrate potential without internal degrees of freedom 



Usir^l,z) = Uoiz) + U,{z) J2 cos(G ■ r,,), (4) 

{G} 

where ry = {x,y) are the coordinates of the Xe atom parallel to the substrate, and {G} is 
the set of the six shortest reciprocal lattice vectors of the substrate. The first term in Eq. 
(1) describes the mean interaction of the atoms with the substrate, and the second term 
describes the periodic corrugation potential. 



Expressions for Uo{z) and Ui{z) were derived by Steele [|T5l assuming that the substrate 
potential Us{r) is a sum of Lennard- Jones potentials between one film atom and all of the 
atoms in the substrate. However, summing Lennard- Jones potentials for Us{r) is not a 
correct description of a metallic surface interaction with a noble gas atom. The corrugation 
potential is reduced (from the value found by summing Lennard- Jones potentials) due to 
electronic screening. For this reason we employ a weaker corrugation potential, as did 



Cieplak et al. in Ref. |jTD|. The corrugation potential we use is 

\ 2" 
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Ui[z) = ae 
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where a = Anexe/AgA*^ / Vs, z* = z / a, a = 2.892 A is the lattice constant of the substrate, 
A* = axe/Ag/o., Qi = 2/^/3. We calculate the Lennard- Jones parameters axe/Ag and Sxe/Ag 
by fitting (i) the position of the minimum of Uq{z) to the distance between a Xe atom in the 
first layer and the ion cores of the substrate {zq = 3.6 A , from |jl6|]), and (ii) the attractive 
well depth to the binding energy of one Xe atom to the Ag substrate {Uo{zq) = —211meV, 
from 13). We find axe/Ag = 4.463 A and sxe/Ag = 13.88 meV. 

Employing Eq. for the substrate potential yields a corrugation amplitude of approxi- 
mately 3 meV. This value is smaller than the corrugation amplitude found using the Steele's 
potential (17.4 meV) [|T^. This last corrugation amplitude leads to slip times two orders of 
magnitude smaller than those found in experiments. We have been able to reproduce the 
results for Kr/Au(lll) [^p!0| employing the above approach. 

Our simulations are carried out at an equilibrium temperature of T = 77.4 °K, and the 
particles move in a three dimensional box of size 20 ax 10a\/3 x 10 a. The time scale for 
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vibrations of the adsorbed film atoms is to = \J {jna"^ /e) = 3.345 ps , with m = 130.1 g/mol. 
Periodic boundary conditions in the x and y directions are employed along with a hard wall 
boundary condition in the z direction at the top of the box. 

We change the coverage by changing the number of Xe atoms A^. We use 60 < < 370. 
All atoms are initially in the gas phase. The atoms condensed in 250 to or less, forming a 
triangular lattice incommensurate with the substrate fcc(lll) surface. A thermostat is used 
only to establish thermal equilibrium. To calculate the slip time we use two methods: 

Method I: An external force (parallel to the substrate surface) is applied to all film 
particles for approximately 100 to. The external force then is removed. This induces an 
initial center of mass velocity Vq in the film which decays at later times as Vo e"*/^. Indeed 
we find an exponential decay as shown in Fig. |l|. The slip-time versus coverage results 
obtained in this manner are shown in Fig. |^. Because the thermostat was turned off while 
the decay of the center of mass velocity took place, the temperature rose by at most 13 °K 
during this process, which occurred for the largest value for the initial velocities used in 
these calculations (0.6cr/to)- 

Method 11: In this method no external force is applied at any time. In experiments using 
QCM 1^, one conventionally describes the data by the linear response in the force per unit 
area 5f to an applied substrate velocity 5V at complex frequency This response defines 
the acoustic impedance of the film, Z{Q = \im.sf^Q{5f / 5V) . The Kubo formula for the 
acoustic impedance has been derived [|12| to be 

Z{0 = -KP2 (l + < \^ C{t)e<'dt\^ , (6) 

where p2 = Nm/A is the film mass per unit area A, and the thermal equilibrium autocor- 
relation function 

Cit - t') = < > (7) 

In Eq. (|^, Vx{t) is the center of mass velocity of the film in the x direction The slip 
time model [|r^ may be written as a Drude-Darcy law impedance Z{Q = —i{(p2)/{l — i(T), 
or equivalently, as the autocorrelation function 



C(t-t') =exp (^-^-^j . (8) 

Equation (^) is the basis of the second method for calculating r. The thermal average 
in Eq. (|^ is calculated performing an average over a time ttot 

< V,{t)V,{t') >= V, {t + s) V, {f + s) ds. (9) 

ttot JO 

In Fig. ^, the values of the slip time as a function of film coverage are shown for the two 
methods described above. The two methods yield results which are qualitatively similar to 
each other and to previous experiments. 

We find a minimum followed by a sharp rise in the slip time for coverages near one 
monolayer. The minimum in r corresponds to the uncompressed monolayer (see Fig. 
We are able to observe directly in the simulations a compression of the monolayer when the 
coverage varies from 0.0563 atoms/A^ to 0.0594 atoms/A^. The interparticle average spacing 
at these coverages are respectively: 4.53 A and 4.4 A. These values are very close to those 
reported experimentally by Unguris et al. ||19[ ( 4.52 A and 4.39 A for the uncompressed 
and fully compressed monolayer). 

The structure factor 

1 ^ 

^(Q) = ]^<Ecos(Q-(r,-r,))> (10) 

has also been calculated. In Fig. ^ we show the Bragg peaks in the Q = {Q^, Qy, 0) plane for 
different coverages: submonolayer, slightly below uncompressed monolayer and compressed 
monolayer coverages. The Bragg peaks for the larger slip times (compressed monolayer) 
are sharper than the Bragg peaks for the smaller slip times (submonolayer and slightly 
uncompressed monolayer). Thus, the lattices corresponding to small slip-times have more 
disorder than lattices with large slip-times. 



Persson and Nitzan |TT| have recently calculated friction for the present system using 
molecular dynamics in a model that includes a Langevin thermostat, which they identify 
with the effect of electronic excitations in the substrate. This thermostat uses unequal 



damping constants in directions parallel and perpendicular to the sliding direction (i.e. 
?7|j = 6.2 X lO^sec^^ and r/^ = 2.5 x lO^^sec^^.) They find a resulting friction constant 
Vtot ~ 6.32 X lO'^sec"^ in their simulations, for a monolayer coverage. Assuming rjtot = 
V\\ + Vphon, one finds that the phonon contribution to friction rjphon/vtotai = 0.02. Thus 
there appears to be a significant reduction in rjphon due to the presence of the electronic 
damping r^y. We have performed simulations using the substrate potential of Eqs. ([4]- 
[5]) with the same Langevin thermostat and damping constants as Persson and Nitzan. 
We find, for a monolayer coverage (0.0563 atoms/A^ ), that Vtotai = 2.45 X 10 sec and 
Vphon = 1-83 X lO^sec"^, so that rjphon/vtotai = 0.75. Thus, although the suppression of the 
phonon contribution to friction is significant, it is not as drastic as that reported in ref. [jlT 



In conclusion, we have performed molecular dynamics simulations of Xe/Ag(lll) in a 
thermostat-free environment. We have been able to qualitatively and semi-quantitatively 
reproduce the experimental data by means of two independent methods which simu- 
lated sliding speeds, which are at least three orders of magnitude apart. In particular, our 
simulation reproduces the distinctive dip in the submonolayer regime found in experiment. 
This result, combined with the fact that the electrical surface resistivity was found to be 
nearly constant in this regime [0 suggests that the phonon contribution dominates over the 
electron contribution for this system. The fact that the simulations yielded comparable slip 
times for greatly varying sliding speeds demonstrate the linearity in velocity of the friction 
governing this system. 
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for his aid. One of us (M. S. T.) would like to thank H. A. Makse for useful discussions. 



7 



REFERENCES 

[1] F. P. Bowden and D. Tabor, Friction and Lubrication (Methuen, London, 1967); J. F. 
Archard, Proc. R. Soc. London Series. A 243, 190 (1957). 

[2] D. Dowson, History of Tribology (Longman, London, 1979). 

[3] J. B. Sokoloff, Phys. Rev. B 42, 760,6745(E) (1990); Phys. Rev. B47,6106 (1993); Wear 
167, 59 (1993); Phys. Rev. Lett. 71, 3450 (1993); Phys. Rev. 51, 15573 (1995). 

[4] J. B. Sokoloff, Phys. Rev. B 52 5318,(1995). 

[5] E. T. Watts, J. Krim, and A. Widom, Phys. Rev. B 41, 3466 (1990); A. Widom and J. 
Krim , Phys. Rev. E 49, 4154 (1994). 

[6] C. Daly and J. Krim, Phys. Rev. Lett. 76, 803 (1996). 

[7] A. Dayo and J. Krim, Proc. 8*'* International Conference on Vibrations at Surfaces, 
Birmingham, England, June, 1996, Surf. Sci. in press. 

[8] J. Krim, D. H. Sohna, and R. Chiarcllo, Phys. Rev. Lett. 66, 181 (1991); J. Krim, 
Comments Cond. Matt. Phys. 17, 263 (1995). 

[9] B. N. J. Persson and A. 1. Volokitin, J. Chem. Phys. 103, 8679 (1995). 

[10] M. Cieplak, E. D. Smith, and M. O. Robbins, Science 265 1209 (1994); Phys. Rev. B 
54, 8252 (1996). 

[11] B. N. J. Persson and A. Nitzan, Surf. Sci. 367, 261 (1996). 

[12] J. Krim and A. Widom, Phys. Rev. B 38, 12184 (1988). 

[13] R. Kubo, M. Toda, N. Hashitsume Statistical Physics II: Nonequilibrium Statistical 
Mechanics, sec. ed. (Springer- Verlag, Heidelberg, 1995). 

[14] J. G. Dash, Films of Solid Surfaces (Academic Press, London, 1975). 

[15] W. Steele, Surf. Sci. 36, 317 (1973). 

8 



[16] G. Vidali, G. Ihm, H.Y. Kim, and M. Cole, Surf. Sci. Reports 12, 133, (1991). 

[17] In our simulations we take Vx as the center of mass velocity minus its mean < Vx >. 
Subtracting the mean center of mass velocity is a standard technique that, in a molecular 
dynamics simulation, ensures that the correlation goes to zero at long times |TB|. 

[18] J. M. Haile, Molecular Dynamics Simulation: Elementary Methods (Wiley, New York, 
1992). 

[19] J. Unguris, W. L. Bruch, E. R. Moog, and M. B. Webb, Surf. Sci. 87, 415 (1979). 



9 



FIGURES 

FIG. 1. Typical decay of the center of mass velocity (in units of cr/to) after applying the external 
force. After the force and thermostat are turned off, the velocity decays; and this portion of the 
graph is fitted to an exponential curve as predicted by the linear friction force law to obtain the 
slip time (method I). 

FIG. 2. Slip time r (in units of sec.) versus coverage (in units of number of particles per A^) 
for the two methods used, (a) Method I. (b) Method II. Solid lines in both figures are three of the 
experimental curves from Ref. |^], included for comparison. There is qualitative and quantitative 
agreement with experiment. Both methods reach a minimum near 0.0563 atoms/A^. 

FIG. 3. Bragg peaks for (a) a submonolayer coverage, 0.040 atoms/A^, (b) slightly below 
uncompressed monolayer coverage, 0.0538 atoms/A^, and (c) the compressed monolayer coverage, 
0.0594 atoms/12. 
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